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ABSTRACT 

o 

An orbiting object in a gas rich environment creates a gravitational density wake 
containing information about the object and its orbit. Using linear perturbation theory, 
we analyze the observable properties of the gravitational wake due to the object cir- 
cularly moving in a static homogeneous gaseous medium, in order to derive the Bondi 
accretion radius re, the orbital distance r p , and the Mach number M p of the object. 

Ph \ Supersonic motion, producing a wake of spiral-onion shell structure, exhibits a single- 

armed Archimedes spiral and two-centered circular arcs with respect to the line of sight. 

'-£5 ■ The pitch angle, arm width, and spacing of the spiral pattern are entirely determined by 

the orbital distance r p and Mach number A4 P of the object. The arm-interarm density 

contrast is proportional to r#, decreasing as a function of distance with a power index 
-*— > , 

of —1. The background density distribution is globally changed from initially uniform 
C^ ■ 

to centrally concentrated. The vertical structure of the wake is manifested as circular 

arcs with the center at the object location. The angular extent of the arcs is deter- 
mined by the Mach number A4 P of the object motion. Diagnostic probes of nonlinear 
wakes such as a detached bow shock, an absence of the definite inner arm boundary, the 



^t 



ON 



Q\ ■ presence of turbulent low density eddies, and elongated shapes of arcs are explained in 

the extension of the linear analysis. The density enhancement at the center is always 
tb/t p independent of the nonlinearity, suggesting that massive objects can substantially 
modify the background distribution. 
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1. INTRODUCTION 

The gravitational interaction of orbiting objects with a background medium is fundamental in 
many astrophysical problems. Objects in gas-rich environments such as planets in protoplanetary 
disks, young stars in natal molecular clouds, substellar mass objects in stellar winds, compact 
objects orbiting around supermassive black holes, and supermassive black hole binaries in colliding 
galaxies are subject to the gravitational interaction with background gas. The orbiting object 



1 Academia Sinica Institute of Astronomy and Astrophysics, P.O. Box 23-141, Taipei 10617, Taiwan; 
hkim@asiaa.sinica.edu. tw 



naturally generates a gravitational density wake in its environment, inevitably followed by the 
reaction of the induced wake as dynamical friction. 

A recent series of studies on the gaseous dynamical friction is devoted to understand the orbital 
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evolution of the object in rectilinear mo tion (|Ostrikerlll999l : ISanchez-Salcedo &; Brandenburg 



2008 



1999 



and in circ ular motion (jSanchez-Salcedo &: Brandenburg 
Kimll2010l ). These theoretical and numerical works have 



covered the effects of binarity, nonlinearity, acceleration, and background stratification. Albeit all 
these factors give rise to appreciable effects, the density wake and the dynamical friction force are 
primarily characterized by the accretio n al an d or bital properti e s of t he object describ ed by the 
linear perturbation analyses in lOstrikerl (|1999l ) and lKim fc Kiml (|2007l . hereafter IKK07I ) . 



Given the importance of the gravitational drag, these linear analyses have been applied to vari- 

ous a s trophysical aspects such as the orbital decay of planets around giant stars ([Schroder fc Connon Smith 



2008 



Villaver &; Livid I2009T) . the prim ordial mass segregation in massive star for ming regions 



(G arav et all2010l ; IChavarrfa et al.ll2010l ). the binary hardening in parent gaseous disk (jBaruteau et al. 
20 111 ) , the deceleration of galaxies as the heating sou r ces of intracluster mediu m against the cooling 



/mil ), tne deceleration or galaxies as the heating sou r ces or intracluster mediu m against the cooling 
flows (JEl-Zant et alJboO^Ecim et al.ll20oi : lKimll2007l : IConrov & Ostrikerlhoosh. a nd the coalescence 



of supermassive black holes in gas-rich mergers of galaxies (jEscala et al. 



2004 



2005 



Dotti et al. 



2006, 



(Barausse 



37) 



2007: iMaver et al J 120071 ). These linear analyses have been extended to the relativistic case 



2007) for extreme mass-rat io inspirals in the vicinity of supermassive black holes. In 



particular, iBarausse fe Rezzollal (|2008l ) found an important difference in the inclination angle of 
the orbits between the dynamical friction and the radiation reaction possibly detectable by Laser 
Interferometer Space Antenna (LISA). 

Also important are the morphological and kinematic features of the density wake per se, as 
they can be directly observed by means of X-rays i n shocked regions, molecu lar line emissions, dust 
continuum emissions, and in dust scattered light. iFurlanetto fc Loebl (12002 ) raised the possibility 
that the nature of dark matter can be examined by detailed X-ray observations of wakes due to 
the motion of galaxies in the intracluster medium, exhibiting distinguishable pict ures between the 
stand ard collisionless cold dark matter and the collisional (fluid) dark matter. ISanchez-Salcedo 
(J2009l ) furthermore proposed that the morphology of a galactic wake in X-ray emission differs in 
purely baryonic Modified Newtonian Dynamics (MOND) from the case of Newtonian gravity with 
dark matter. T h ese au thors only inspected the local attributes of the gravitational density wakes 
using lOstrikerr s (|1999l ) formulae for the cases of the perturbing object in rectilinear motion. 



However, many observable objects have fairly short orbital periods compared to the dis persion 
timescale of the background gaseous matter, leading to global shaping of the induced wakes. LKK07 
showed that the curvilinear orbit modifies not only the curvature of the wake but transformes the 
overall structure entirely because of the memory effect^. One application of the large scale wakes 
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is to a case that planets or brown dwarfs lie in cool giant envelopes, whose wakes are possibly 
observable at large distances from their bright parent stars. Future high sensitivity/resolution 
observations of these wakes can provide evid ence of the unseen substellar mass objects, further 



see 
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all owing us to derive t heir orbital properties (|Kim &; Taamll201ll in prep.; for a brief introduction 



Although IKK07I calculated the density wake, they primarily focused on the dynamical friction 
force that affects the initially circular orbit and the description of the general wake feature was 
somewhat superficial. A more detailed understanding of the wake characteristics is necessary 
in order to provide an interpretative framework for both further theoretical works under more 
complicated situations and future observations to de tect th e gravitational wakes of hidden objects. 



In this paper, we extend the semianalytic work of IKK07I highlighting the observable quantities, 
and determine the detailed structure and kinematics of the gravitational wake due to an object in 
circular motion. Particularly, we study the density structure in great detail based on an analytic 
linear analysis, which also helps to understand aspects of the simulated density wake in a nonlinear 
regime. We provide the density contrast of the wake as a function of the object mass, orbital 
radius, orbital speed, and the background sound speed, providing a way to obtain the information 
of unseen orbiting object from the observation of the wake at a large distance. 

This paper is organized as follows. In $2] we review the derivation and findings in previous 
works as the starting point of this investigation. In Sj3] we present the nature of linear wake as 
revealed in the shape ( §3.11) . the density enhancement ( §3.2|) . and the velocity distribution ( §3.3[) . 
provided that the weak gravitational influence is thoroughly described in the linear perturbation 
analysis. In ^21 based on the linear wake properties, we suggest a new explanation for nonlinear 
effects manifested on the wake when the perturbation in discontinuity is saturated. In £}5] we 
summarize the newly found properties of the gravitationally induced wake of a circularly orbiting 
object in a static uniform gaseous medium, and briefly consider applications to a few astronomical 
objects. 



2. BACKGROUND REVIEW 

Starting from the basic equations of hydrodynamics with an external time-dependent gravita- 
tional potential & p but excluding self-gravity of the gaseous medium, one can obtain a linearized 
three dimensional wave equation for the density enhancement a = 5p/p under the assumption of a 
small perturbation (a <C 1), 

^2 1 9 2 a 1 _ 2 _ ,„. 

^-SJF-a 7 ^ (1) 



of losing their maximum energy capacity after being repeatedly recharged to partially discharged states. We apply 
this term here for the wake repetitively perturbed by the revolving object over the former perturbations remained, 
resulting in a loss of density contrast, as we will mention at the end of this paper. 



Here, c s = ("fp/p) 1 ' 2 is the so und speed of un perturbed gas with density p, thermal pressure p, and 
adiabatic index 7 in general. lOstrikerl (|1999l ) applied it to the situation that the perturbing object 
is in rectilinear motion thr ough an infinit e homogeneous gaseous medium, utilizing the retarded 
Gree n's function technique (| Jackson! Il975l ) based on the Fourier transform (see also I Just &: Kegel 
1990). A purely analytic solution for the density enhancement was derived, which can be expressed 
as 

I>(^t-d|[) (2) 
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in the axisymmetric coordinates (d\\, d±) with the origin at the perturbing object, d\\ along the 
tail (i.e., in the opposite direction of the orbital motion) and d± in the orthogonal direction. Here, 
we define the most important parameters in this paper, rs and M p , as the Bondi accretion radius 
and the orbital Mach number of the object, respectively. In the summation of the Heaviside step 
functions H, the primed distance d!, denotes the original position of the gravitational perturbation, 
initiated by the object and currently arriving at the observing point (dn, d±_). Thus, the integrand 
means that the perturbation sources should be on the trace of the object, inevitably not able to 
be further than the initial position of the object, V p t > dl. The summation of only the step 
functions indicates the exactly same contributions from the individual perturbations gathered at 
a time. Accordingly, it can be substituted simply by the number of perturbations at the point, 
N = Yld' H (^p^ ~~ d!i), which divides the space into three regions showing discontinuous density 
and velocity distributions: AT = 2, 1, and in the cone region, in the ice cream region, and outside 
of these two regions (as following the author's original nomenclature). These three regions are 
separated by the sharp edge of a Mach cone and the modest edge of a sonic sphere about the initial 
position of the object. As the orbital Mach number Ad p decreases, the apex of the Mach cone 
attached to the object approaches to, and finally disappears in, the sonic sphere with the transition 
occurring at the sonic speed (M p = 1), eliminating the shock front from the density wake. 



Using the same retarded Green's function technique, IKK07I applied the three dimensional wave 

equation (eq. [T|) to the linear wake induced by a circularly orbiting object at distance r p from the 

system center in a static uniform gaseous background. The resulting density enhancement is given 

as 

t_b_ v-^ H (V p t/r p - (p n 
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in cylindrical coordinated (r, (p, z) with the equatorial plane aligned to the orbital plane of the 
perturbing object. All the possible perturbations must be launched from the object passing through 
the locations (r p , ip' , 0), where <p' satisfies the condition, 

1/2 
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2 ip is the angular distance from the object in the opposite direction of the orbital motion, tracing the induced 
w ake tail , i.e., the perturbing object has moved from V p t/r p (>0) to 0. This is corresponding to — s in the notation 



of KK07 



This relation describes that the individual perturbation propagates from the launching point 
(r p , ip' , 0) to the observing point (r, <p, z) as a sound wave during the time that the perturb- 
ing object travels the angular distance r p ip' along the orbit with the speed Ai p . In the case of 
curvilinear motion, the summation is not only restricted to the step function, implying different 
contributions from the different elements of the perturbations. This diversity of contributions is 
the primary reason of the complications arising in a density wake of the circularly orbiting object. 



KK07I accurately calculated the wake density in equation ([3]) by computing the positions of 
perturbation origins in the nonlinear equation §£§ using root-finding algorithms based on either 
Newton-Raphson or bisection iterations. Their accurate numerical values, however, do not provide 
deep physical insight for the global features of the induced density wake , which are useful to probe 



the orbital properties of the perturbing object. In this paper, we revisit iKKOTT s work on the linear 
density wake of a circularly orbiting object by modifying their analytic formulae and by inspecting 
the global distributions of the calculated density enhancement and velocity fields, if necessary, in 
comparison with the linear trajectory counterparts. In addition, we have performed high resolution 
three dimensional hydro dynamical simulations using an adaptive mesh refinement code, FLASH 



(JFrvxell et al.ll2000l ) , to investigate the response of a gaseous medium to the influence of an extended 
object, characterized by a gravitational softening radius r s in the range of 0.01-0.1 r p (alternatively, 
0.1-25 re) for Ai p = 0.2-10. This parameter space includes the nonlinear regime showing diff erent 
behaviors of the wake from the ones in the linear analysis (see also lKim fc: Kimll2009l ; lKimll2010l ). All 



the simulations treat only a pseudo-isothermal gas with 7 = 1.00001 in the FLASH code, providing 
a constant sound speed c s in the entire computational domain. 



3. LINEAR WAKE 

3.1. Morphology 

An object orbiting with a subsonic speed (A4 P < 1) creates a smoothly distributed wake 
without a shock, showing a gradual decrease of density along the distance from the perturbing 
object. The isodensity surfaces of this subsonic wake have the shapes of oblate spheroids with 



ellipticity e = M p near the perturbing object (see Appendix A of lKK07l for derivation), gradually 
twisted in a trailing fashion with the low-density outermost shells tending toward comma shapes. 
The morphology of a subsonic wake is simple and easily understood since the wake receives only one 
sonic perturbation everywhere, contrary to the supersonic wakes (M. p > 1) investigated in depth 
in this section. 

The wake induced by a circularly orbiting object with supersonic speed has a spiral-onion shell 
structure in three dimensional space, having distinct boundaries for the high density region and 
the lower density inter-structure regions. The density contrast is primarily caused by the piling 
up of the delayed gravitational perturbations inside the high density structure. Equation §S§ gives 
the positions of the origins, (p', for the perturbations reaching the observing point (r, ip, z) at time 



t. After some algebra from equation dH), IKK07I showed that the high density structure possesses 
N = In + 1 perturbations, where n is an integer number depending on the orbital Mach number 
such as 1 (M = 3) at M p ~ 1 - 4.6, 2 (M = 5) at M p ~ 4.6 - 7.8, 3 (A/" = 7) at A* p ~ 7.8 - 10.9, 



and so on (see Appendix B of IKK071 for a more detailed discussion). For comparison, in subsonic 
cases {M p < 1) the wake receives one perturbation (A/ - = 1) everywhere, corresponding to n = 0. 



Extending the analysis in Appendix B of lKK07l . in which the number of perturbations of the 
arm structure in the orbital plane is obtained, we derive a purely analytic formula to describe 
the exact shape of the enhanced density structure in three dimensional space. The details of the 
derivation are omitted but the essence is to use the property at the boundaries, where the number 
of perturbations changes so that the curve on the right hand side of equation (jj]) has the same slope 
as the line on the left hand side. We derive an expression for the shape of the enhanced density 
structure given by 
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for the outer (upper sign) and inner (lower sign) boundaries, where X and Y are defined as X 
A4 p (d++d~)/2 and Y = A4 p (d+—d~)/2. Here, the normalized distances d— = ((r 



+ z< 



,1/2 
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from the perturbing object and d+ == ((r + r p ) 2 + z 2 ) jr v from the mirror point about the sys- 
tem center. The solutions of these formulae determine the boundaries of the high density region, 
revealed in Figure [T^i-b as a single-armed spiral in the orbital plane (z = 0) and arcs in a meridional 
plane (tp = 0). 



3.1.1. Spiral 

Shaping a spiral arm having a constant width is physically not so trivial for the density wake 
in a uniform background. Because of the gas pressure, the gaseous matter accumulated behind the 
object tends to disperse with the sound speed in a static uniform background. Thus, the excess 
density region should form an extending tail similarly to the Mach cone shape in the linear trajectory 
counterpart, unless the background density configuration is changed. As we will show in §3.2.11 the 
initially uniform background transforms the distribution into a hydrodynamic equilibrium about 
the system center. As the result of the centrally concentrated background, the inner part of the 
spiral does not propagate inward but outward, following its outer part. 

In this section, we construct an explicit form of the spiral shape, which is important in order 
to derive the pitch angle, the arm width, and the spacing of the spiral pattern, relevant to the 
object orbital parameters. Equation © at z = yiel ds the p olar equations for the boundaries of 
the spiral arm in the orbital plane (see also eq. [B2] in|KK07]), 
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at z = 0, (6) 



in which the upper and lower signs are for the outer and inner bo u ndari es, respectively. The outer 
boundary of the arm, originated from the Mach cone in lOstrikerl (119991) , lies not only outside the 
orbit (r/r p > 1) but it also extends to r/r p = Ai- 1 (see Fig. 1 in IKK07I ) . At this innermost point 
of the outer boundary, the sonic sphere due to the influence of the gravitational perturbation from 
the object at the initial location starts to overlap with the curved Mach cone tailing the object, 
forming the inner boundary of the spiral arm. 

In order to see the global features of the wake in the orbital plane, we differentiate the com- 
plicated forms (eq. [6]) following the exact shapes. Using the angle transformation formulae for 

cosine, it is reduced to 
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for both outer and inner arm boundaries in the orbital plane, except for a small part of the outer 
boundary inside the orbit (.M" 1 < r/r p < 1) having the opposite sign. We define an angle as its 
tangent value, tan0, corresponding to equation ([7]), which leads to the alternative form 
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resulting in the familiar result = sin 1 M. 1 at the object p osition r/r p = 1 as the half-opening 



angle of the curved Mach cone (JKK07I ; see also lOstrikerl Il999l ) . This angle decreases along the 



distance r/r p as seen in equations ([7|) and flSJ), but it converges very quickly to G = tan" 1 .M" 1 , 
showing a small deviation of < 1° at r/r p > 2 for A4 P > 2 (r/r p > 4 for Ad p = 1). This quasi- 
constancy of suggests that the density wake in a static gaseous medium indeed approximates an 
Archimedes spiral, satisfying r/r p = (pM~ l + C where the constant C differs between the inner 
and outer boundaries of the high density arm region. Figure [Ik. displays the density enhancement 
in the orbital plane (z = 0) of the perturbing object, and exhibits a high density region enclosed 
by two Archimedes spirals shown in dashed (C = 1) and dotted (C = —0.7 for A4 P = 10) lines. 
The constant C's are calculated at r/r p = 1 using equation ([6]); the value for the inner boundary 
changes with A4 P , while it is unity for the outer boundary. 



3.1.2. Arcs 



An investigation of the wake shape in a meridional plane is needed to provide an interpretative 
framework for edge-on observations. Compared to the asymmetric spiral pattern seen in the face- 
on line of sight, this arc pattern is not properly spotlighted and hence not analyzed well. In this 
section, we show the purely analytic solution for the arc shape of the gravitational wake, and use 
a geometrical approach to fully understand the wake shape relevant to the object properties. 

A trivial solution of equation ([5]) is both X = 1 and Y = 1 in the meridional plane wherein the 
perturbing object lies ((p = 0), which indicates the (vertically) highest intersecting points between 



the outer and inner boundaries, i.e., the extension limit of the intensively perturbed region. X = 1 
and Y = 1 stand for an ellipse and a hyperbola, by definition, with the common foci at the location 
of the perturber and the mirror point about the orbital axis. However, both X = 1 and Y = 1 
actually correspond to the same formula for a hyperbola 

M 2 
M 2 p (r/r p ) 2 -j^(z/r p )* = l, (9) 

since we only consider the shape of a wake induced by a supersonic perturber (M p > 1). The 
vertex of the hyperbola, r/r p = Mp 1 , corresponds to the distance of the junction between the 
outer and inner arm boundaries in the orbital plane. The extension limit of the high density arcs 
in the vertical direction is simplified to the asymptotes of the hyperbola, 

z = ± {M 2 p - 1) 1/2 r, (10) 

which are overlaid by black lines in Figure [lb. This indicates that the angular size of the arcs 
is 2 tan -1 (.Mp — l) 1 ' 2 , further extended with higher orbital Mach number. This vertical stretch, 
however, does not depend very strongly on M p except for near unity (i.e., 60°- 90° for M p > 2). 

In Figure [lb individual segments of the boundaries, satisfying equation ([5]), appear to be 
circular arcs about the centers at either the object position or the mirror point. In the plane of the 
object (f = 0; x > 0), the arcs extended from the outer boundaries of the spiral arm in Figure [TK 
(i.e., x/r p = 1, 1 + 2-jrMp 1 , 1 + 4irM p \ 1 + 6irM p l , 1 + 8irM p \ and 1 + WirM p l at z/r p = 0) 
center around the perturbing object at (x,z)/r p = (1,0). These concentric arcs have endpoints 
on the hyperbola defined in equation Q (or, the asymptotes in eq. [TO]), except for the point 
(not visible) at the perturber position and the curve passing through (x,z)/r p = (1 + 2irMp 1 ,0), 
which completes its full circle. The inner x-intercept of this circle at x/r p = 1 — 2ttMZ 1 actually 
corresponds to the inner edge of the curved Mach cone, which is considered as a part of the outer 
arm boundary (see §3.1.1[) and, here, this is clarified. On the other hand, the arcs crossing the 
inner arm boundaries on this side (x > 0) have the center at the mirror point (x,z)/r p = (—1,0), 
resulting in larger radii compared to the radii of the outer counterparts. The outer boundaries on 
the opposite side (x < 0), i.e., the arcs having intercepts at x/r p = — 1 — irMp 1 , —1 — 37rM p , 
—1 — 5ttM p 1 , —1 — 7irM p 1 , and —1 — 97rM p 1 , are concentric about the mirror point of the object, 
(x, z)/r p = (—1, 0). On this side (x < 0), instead, the center of the inner boundaries is at the object 
position, (x,z)/r p = (1,0). 

The inner and outer boundaries, in fact, constitute the same circles, which is better visualized 
in Figure [2] for the number of perturbations. Particularly, the largest sphere about (x,y,z)/r p = 
(1, 0, 0) limits the region already in a steady state, receiving 1, 3, 5, and 7 perturbations in different 
localities. The arcs (or circles) in Figure [2j) are formalized by equation ([4|o with the sources at 
(p' of 7r-multiples in the plane of ip = 0. The perturbation launched every time the object returns 



Eq. (0) gives the same results. Since \(X 2 — l) 1 ' 2 =p (Y 2 — l) 1 | in this equation is originated from tp' in eq. Q, 



back to the original position (p' being a multiple of 2n) propagates as a sound wave leaving the 
vestige on the spherical shell with the radius of multiple of 2irA4~ ; the perturbation generated 
half an orbit (it in angular phase) later is marked on the sphere, smaller by txM~ x , with the center 
at the mirror point in the same meridional plane. And by extension, the spiral in the orbital plane 
(Fig. [2k.) emerges as an aftermath of the phase difference of the perturbation source. 



3.2. Density Enhancement 

The density enhancement, a, in the pattern follows a complicated rule because of the ac- 
cumulation of gravitational perturbations launched from the object at different positions on the 
curvilinear orbit. This complication is particularly caused by the different contributions of the re- 
spective perturbations, unlike the case when the object is in a straight-line trajectory. The density 
enhancement profiles are displayed in Figure [Tfc as a function of radial distance in various directions 
along the spiral arm (the lines in red, yellow, green, and blue colors) and along the orbital axis 
(the cyan line), in which all the curves appear well located within the black lines (see below). The 
profiles in the orbital plane reveal the highest density peaks at the object orbit, r/r p = 1, and the 
peak values decrease with distance from the orbit. The minimum density enhancement in azimuthal 
directions, i.e., the baseline of the group of profiles, gradually decreases following the bottom dashed 
line (up to r/r p = 0.4) starting from a = rs/rp at the orbital center. At r/r p = 0.4, however, the 
baseline jumps to, and follows, another dashed line, and another jump occurs at r/r p = 0.75. In 
this example case of Ad p = 10, the baseline moves its paths among the dashed lines through three 
steps before reaching the orbital radius. Beyond the orbital radius, it follows the third dashed line 
from bottom until it abruptly drops at r/r p ~ 2.5 apart from the path. This last drop of minimum 
density enhancement coincides with the sudden disappearance of the inner arm boundary, which 
is related to the orbital evolution time of the object. Using the spiral functional forms deduced in 
§3.1.1} it is shown that the outer and inner arm boundaries reach r/r p = 4.1 and 2.4, respectively, 
at the moment that the perturber completes 5 orbits (i.e., <p = 107r). The inner arm boundary 
lags behind by 2.7-turn since it starts from (r/r p , p>) = (A4" 1 , 1 + 0.7A^ P ), deeper inside than the 
starting point of the outer boundary at (r/r p , p) = (1, 0). 



eq. ([5]) is equivalent to (X 2 + Y 2 — 2XY cos(<p' — if)) 1 ' 2 = f' ■ I n the meridional plane of ip = 0, the gravitational 
perturbations from the object located, in the past, at ip' — 2-7TX0 (current), 1, 2, 3, 4, and 5 (initial) form the arcs 
as the segments of circles, rf_ = (X — Y)A4 p ~ 1 = 2-kM~ v 1 x0, 1, 2, 3, 4, and 5. Similarly, <p' = irxl, 3, 5, 7, and 9 
correspond to the concentric arcs on d + = (X + Y)A4p 1 = irMp Xl, 3, 5, 7, and 9. 
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3.2.1. Interarm Region 

Equations ©-(JU) indicate that the minimum of density enhancement in the azimuthal direction 
satisfied 

r _B nn 

" mln ( r 2 + z 2 + r 2)l/2 + rA1p - ^ 

This is plotted in Figure [Tb by the bottom dashed line, outlining the interarm density enhancement 
for all directions. The incremental steps, observed above, follow the dashed lines exhibiting a = 
amin(^) z = 0) X 1,4, 7, ■ • • , 1 -+- 3ra, where n = (W — l)/2 i s 0, 1, 2, etc. depending on the orbital 



Mach number A4 P (for n see N3.ll originated by eq. [B3] in lKK07l ). 



To examine the dependence of the minimum density enhancement on Ai p (technically, n), we 
compare the a-profiles for M. p of 0.8 (n = 0), 2.0 (n = 1), 4.5 (n = 1), and 5.0 (ra = 2) in Figure O 
having n + 1 dashed lines in each panel. Unlike the subsonic case (n = 0) having smooth a-profiles 
above a m i n in Figure [3h., the supersonic cases (n > 1) in Figure [3b-d show arm structures with 
sharp boundaries of higher density than in the middle of the arm. In Figure [3b for A4 P = 2 (n = 1), 
the interarm density is still represented by a m [ Q (lower dashed) and the local minimum inside the 
arm is higher than 4a m i n (upper dashed). Increasing M. p (up to 4.6) increases the arm width with a 
corresponding decrease in interarm extent (see Fig. [3b). The minimum arm densities concurrently 
approach 4a mrn . At the critical Mach number M p = 4.6, the outer boundary of the arm meets 
with the inner boundary of the arm one turn ahead, leaving no gap between arms. Because of this 
overlap (outside the orbit), the colored lines in Figure [3jd for A4 P beyond 4.6 no longer exist on 
the bottom dashed line (a = a m ; n ). Instead, the minimum density is now given by the value of 
a = 4a m ; n , characterized as the local density minimum within the (infinitely) broad arm, above 
which a new narrow arm forms. The new arm possesses an increased number M of perturbations, 
changing n from 1 to 2, because this arm is basically the overlapping region. Increasing A4 P above 
4.6 repeats the same broadening of the new arm until M p reaches another critical Mach number, 
7.8. Above this critical Mach number, an overlapping event causes another jump of the minimum 
density to (1 + 3(n — l))a m in with the increased n (>1, in general), and leads to the density of the 
newly formed arm above (1 + 3n)a m i n . 

As shown in Figure [lb, the density structure at high latitudes (\z\ > 0) is more complicated 
because of the partial overlaps of overdense areas closed by arcs with different radii. The overdense 
regions induced by a slower object (1 < A4 P < 4.6) are composed of simple crescent-shapes without 
overlap. However, once overlapped (A4 P > 4.6), the regions of interest exhibit enhanced densities 
because they swallow the perturbations of each crescent-shaped structures. For instance, the case 



4 Strictly, eqs. ([3]|-(|4} imply the minimum value of density enhancement in the central part r/r p < 2M p 2 is 
better represented by a ~ tb / ((r + r p ) 2 + z 2 ) than « m i n defined in eq. {TT}. This is simply seen in an extreme 
case of M v = so that 2M~ goes to infinity; its density wake in hydrostatic equilibrium about the static perturber 
at (r, if, z) — (r p , 0, 0) shows the density enhancement a — exp(rs/d) — 1 fs rs/d (when r^/d ^C 2) at the distance 
d = (r 2 + z 2 +r 2 — 2rr p cos ip) 1 ' 2 from the perturber, which has the minimum value at (p — n. 
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in Figure Q] (A4 P = 10) experiences 3 steps of density enhancement (corresponding to the number 
of perturbations W = 7) within the discernible spiral arm in the orbital plane. But the number of 
perturbations decreases at high latitude where the overlapping events are reduced (see also Fig. [2b). 
Eventually above the extension limit (Eq. [9]), the gas flow experiences only one perturbation so 
that the density enhancement has the value of a m i n (r, z) (see cyan line in Figure [lb along the dotted 
line, a = a m \ n {r = 0, z)). 



3.2.2. Arm Region 

In order to understand the distribution of the peak densities along distance, it is preferable to 
refer to the simpler case of the perturbing object moving on a linear trajectory, w hose resp ective 



perturbations equally contribute to the density enhancement of the induced wake. lOstrikerr s for- 
mula for the linear motion of a perturber (eq. [2]) has roughly an anticorrelation with respect to the 
distance fr om the p erturbing object with the amplitude of tb l-M 2 , — lj" 1 ' 2 per each contribution. 



Modifying [OstrikerJ's formula, for circular orbit cases we hypothesize a density enhancement due to 



one perturbation, 

Tg 1 

ai= \r-r p \\Ml-l\^ (12) 

in the orbital plane. However, in reality, the circular motion modifies the density enhancement 
from a simple collection of a\ by the number of perturbations. Particularly, gas pressure at the 
central region prevents the propagation of Mach waves, forming instead an equilibrium state with 
the central value a = rs/rp. A rough guess of a\ can be used in an empirical formula for the peak 
density enhancement in the spiral arm pattern, 

a £3 a\ + 2A/a m i n , (13) 

which shows fairly good fits in Figures [lb and [3] (see black solid curve), except at the exact locations 
of the arm boundaries where infinite values are introduced by the point mass perturber (see below). 



3.2.3. Effect of Object Size 

Although, mathematically, the perturbation from a point mass causes infinite density at the 
exact locations of arm boundaries, it is limited to the immediate vicinity of the arm boundaries, 
detectable only at extremely high angular resolution. For an extended object, the very sharp 
density peaks smear out depending on the size of the object, r s , but the peak densities ch ange 



signi ficantly only within a distance of the order of r s from the boundaries (see Fig. 6b of 



Kim 



2010i ). In order to see the effects of the object size on the global structure of the density wake, we 
compare the density enhancement a in Figures [IJ U and [5] for the cases of r s /r p = 0, 0.01, and 0.1, 
respectively, where the object size r s is defined as the gravitational softening radius of a Plummer- 
type object. From Figures [1^,-b, S^i-b, and Ek,-b, one finds that the sharpness at boundaries are 
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reduced with the softened gravitational potential, while there are negligible changes in the global 
shape and morphology of the wake. The profiles of a along the spiral pattern in Figure Ht show 
the smoothing of arm boundaries compared to Figure [lb, chiefly by reducing the peak densities. 
Figure [5b indicates, however, that the minimum value of a is also changed in the case of the further 
extended object with size r s comparable to or larger than half of the interarm spacing for the point 
mass counterpart, reducing the contrast between arm and interarm densities. 



3.3. Velocity 

In Figure [6] zooming in on a high density structure of the wake, the initially static medium 
is revealed to attain an approximately spherical velocity field. In order to understand the velocity 
features of the wakes for circularly o rbitin g objects, we first inspect the linear trajectory coun- 



terparts formulated in iKim &; Kind 's (]2009i ) equations (A19)— (A20). From their formulae for the 
velocity components Vji and V± parallel and perpendicular to the object motion, respectively, it 
is straightforward to show that in the linear Mach cone (half-opening angle of @m = sin -1 A4" 1 ) 
(1) the perturbed gas flows toward the edge of the cone perpendicularly, i.e., the fluid angle is 
tan~ 1 (|Vj_/Vj||) ~ 90° — @m hi th e vicinity of the boundary, deflected f rom 0° headi ng for the 



objec t behind it (see also Fig. 1 in ISanchez-Salcedo &; Brandenburg! Il999i and Fig. 8 in iNamoun 



20 111 ): and (2) the fluid speed in terms of Mach number approaches the value of the density en- 



hancement (M. ~ a) close to the boundaries but decreases to M. = rg (Af p d) * alon g the line 



of the object motion at distance d from the object (see also Fig. 1 in lKim &; Kimll2009l ). Overall, 
the fluid in a Mach cone has the maximum velocity near, and toward, the boundaries, although it 
deflects in angle for considerably slower fluid parts giving smaller contribution on average. 

In the extension of this tendency of the flow to aggregate perpendicularly to the shocked 
boundary, it can be understood that the velocity vectors in Figure [6] appear to be spherical especially 
around the tightly wound spiral- and circular arc-shaped structures. The spiral boundaries in the 
orbital plane have a small pitch anglqj of tan _1 (A4 p r/r p )~ 1 (see t j3.1.ip < 10° at r/r p < 3 for 
A4 P > 2 up to at r/r p < 6 for the extreme case of M. v = 1, making the velocity vectors further 
spherical at larger distances even in the marginally supersonic case. The velocity fields normal to the 
circular arcs in the meridional planes are obviously spherical at large distances, notwithstanding 
that the centers of the arcs are technically not at the system center but at (r,z)/r p = (±1,0) 
( §3.1. 2j) . As implied by the above inspection of the velocity field in a Mach cone, the flow directions 
inside the high density structure are not normal to the boundary shapes. However, the magnitudes 
of speeds are negligible compared to the maximum speeds around the boundaries, KA = re.M~ 1 |r — 
r p \~ l x 2 (r^/rp)" 1 ' 3 in Mach number, which is empirically found in the range of A4 p = 1-10 and 
r s/ r p = 0.01-0.1 with rs/fp = 0.01. It also deserves mentioning that the flow aggregates toward 
the outer boundary, while it tends to diverge from the inner boundary originated from the behavior 



The pitch angle of a spiral is defined by tan 1 (--jr-), which at r — r v corresponds to Q defined in £13. 1.11 
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of flow on a sonic sphere. 



4. NONLINEAR WAKE 

The density enhancement of a wake is related to rB/r p . As this value increases, for a fixed 
M. p , the arm density increases a nd finally becomes saturated, showing unexpected behaviors. In 



previous numerical simulations (JKim &: Kiml 120091 ; iKiml |201fJ) , the nonlinear features such as a 



detached bow shock, turbulence, a disappearance of the inner arm boundary, and an increase of 
background density were exhibited without a clear explanation. In this section, we describe the 
nonlinear features of the wake based on the linear analysis in the previous section, emphasizing a 
new explanation for the detached bow shock as a phase shift of the inner boundary (defined for 
the corresponding linear wake). With this new explanation, we show that the nonlinear features 
are natural consequences, thus understandable in terms of an extension of the properties of linear 
wakes. 

In Figure El we show the vertical structures of wake density and velocity field for M p = 5, 
just above a critical Mach number 4.6, causing the shallow overlaps between the crescent-shaped 
overdense regions predicted in the linear analysis. The wake shape in FigureEH (for rs/rp = 0.01) is 
well explained by the linear prediction in §3.1.21 (dashed and dotted lines) , except for the broadening 
of the edges to form high density belts of width 2r s . In addition, the inner boundaries of the 
crescents (dotted lines) are enlarged by r s in radius compared to the prediction, since they form 
when the sonic sphere encounters the curved Mach cone broadened by r s . The velocity field in 
Figure EH shows the flow from the inner dotted arcs to the outer dashed arcs. These dotted and 
dashed lines are duplicated in Figure Eb as a reference to the flow pattern in the linear wake. 

The wake in Figure Eb is induced by a 100 times stronger gravitational potential than in 
Figure Eh- This nonlinear wake preserves the density discontinuities at the outermost edges of the 
crescent shapes in the linear counterpart (dashed lines), however they no longer exhibit a belt-like 
appearance. The expected inner edges of crescents (dotted lines) are absent in the nonlinear wake, 
which are in fact shifted outwards beyond the outer edges. For instance, the high density band in 
Figure EH near (x,z)/r s = (—1.0, 1.5) (denoted by "A") tends to fill out the low density region at 
smaller z/r p as well as to supply additional material to the outer boundary band at larger z/r p , 
achieved only when the perturbation is sufficiently strong (Fig. Eb)- The thermal gas pressure close 
to the crescent's apex ("A") is sufficient in nonlinear regime (Fig. Eb) to spread the material to the 
point "B" beyond the expected outer boundary. But the flow originated at a lower latitude ("I") 
is unable to reach the outer edge of the crescent because of the large distance. Instead, it turns 
back to the inner part, shedding material on the way ("II"), and eventually reaches the edge of the 
smaller crescent defined in linear regime ( "III" ) . This process raises the density in between the two 
shock fronts, forming a hill shape of density discontinuity (passing through the point "II"). 

An extreme nonlinearity is exhibited in Figure [8] for the case of a slower orbital motion (A4 P = 
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2.2) even with a smaller gravitational impact {rs/fp = 0.5) than in Figure[7b. In fact, the degree of 
nonlinearity can be defined as t\b = rs/r p (A4p — l)" 1 for the wake driven by a circularly orb iting 
object at distance r p , and the criterion r/B ~ 0.1 for distinct nonlinear behaviors is found by 



Kim 



( 2010l rl. In Figure [8^, the high density region in the orbital plane is not enclosed within the 



dashed and dotted lines representing the spiral arm boundaries of the corresponding linear wake. 
Instead, the sharp discontinuity matches the solid line, which is just the inner boundary of the 
linear counterpart (dotted line) shifted in phase by it. This phase shift is mainly caused by the 
high pressure at the junction (r/r p = A4" 1 ) between the expected outer (r/r p = —ipM.~ x + 1) and 
inner (r/r p = (pM^ 1 + 0.2) boundaries, if the nonlinear effect is ignored, forming a sharp tip. The 
matter at the tip tends to spread the high density over a larger region due to the steep pressure 
gradient, unfolding the inner boundary to match the curvature of the outer boundary. This phase 
shift would be saturated to the maximum value of 7T in the very nonlinear case, forming a mirror 
reflecting curve of the inner spiral boundary defined before the nonlinearity is manifested. 

The solid curve in Figure [8^ can be easily mista ken for a p hase shifted outer boundary of 
a linear wake (e.g., 75° for Ad p = 2 case in Fig. 3f of iKiml 120101 ) instead of originating from its 



inner boundary, because the shapes are identical between the two spirals. The 7r-shifted inner 
boundary of the linear wake corresponds to the outer boundary of the linear wake shifted by 
it + 2(M 2 p - 1) 1/2 - cos" 1 (27W- 1 - 1), which gives 80° in advance for M p = 2.2. Although it 
is not clearly seen in the orbital plane, a density discontinuity is found at the location of the 
outer boundary (defined in the linear counterpart) in Figure [8b for the vertical structure, i.e., 
the arc structure having x-intercept at ~ 4r p (see also Fig. [7b). This presence of the expected 
outer counterpart suggests that the origin of the discontinuity along the solid line is not the outer 
boundary but it is likely interpreted as the modified inner boundary of the linear wake with the 
above explanation. As a result, the interior of the nonlinear wake is filled up without a definite 
inner boundary, except for a moderate interference near the black solid lines in Figure [8b (eq. [10]). 
The vertical structure of the nonlinear wake is more elliptically shaped than circularly shaped. 

As the extension of the solid curve in Figure [8^., the material immediately behind the object 
also flows over the curved Mach cone shape, advancing 6° .7 in angl e along the orbit. This angle of 



the isothermal wake (7 = 1) excellently matches the fit in lKiml (j2010l ) for the shock standoff distance 
in adiabatic gas (7 = 5/3). It implies that the shock detachment is a common nature of nonlinear 
flows irrespective of adiabatic index 7 as long as the sonic speed is defined as c s = (jp/p) 1 ' 2 - 
In the process of the overflow creating a new shock front in advance, the flow is highly unstable 
between the old and new shock f ronts, generating low density eddies in extremely nonlinear cases. 



According to iKim fc Kiml (J2009I ) for linear trajectory cases, the low density eddies survive only 
when the detached shock distance from the object in the lateral direction is larger than vbM~ 2 . 



6 For rectilinear motion of the object, the similar nonlinear features are characterized by the definition of the 
degree of nonlinearity with r p replac ed by the gravitat ional size r s of the object, ija = tbJTs (Mp — 1) , and the 



criterion for nonlinear regime is ~ 2 ( Kim fc Kimll2009r ) 
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It is checked that the lateral distance of the shock in Figure [Hh. is also larger than this value, 
corresponding to 0.1r p . 

Lastly, Figure^ shows the density enhancement of the wake in this extremely nonlinear regime. 
The peak values do not have a systematic change from the linear prediction (solid curve), while the 
minima do not match any of the dashed curves unlike in the linear regime. The minima instead 
increase to the dotted curve, a = rs (r + ?"p) _1 , reducing the arm-interarm density contrast in the 
entire wake. The wake due to circular motion attains a hydrodynamic equilibrium state about the 
system center. In contrast, the wake due to rectilinear motion has the equilibrium center about 
the object itself, so that the amplitude of gravitational influence is scaled by rg/r s relevant to the 
object size r s . The wake of a circularly orbiting object actually tends to achieve quasi-equilibrium 
about both the object and its orbital center. However, the former influences the localized area at 
distance from the object of at most rs{Mp — 1) _1 (if this distance is greater than r s , and M p > 1; 



Kim &: Kiml 12003 ; lKimll2010l ). while the latter affects the wake properties globally. The central 



value of density e nhanc e ment maintains the value of tbJt v independent of the nonlinearity, thus 
readily explaining iKimr s (120101 ) claim that the background density along the o rbit i s proportional 
to this value as numerically found a ~ 0.46(rs/r p ) 11 in equation (10) of lKiml ()2010l ). 



SUMMARY AND DISCUSSION 



The formation of a gravitational wake is a characteristic feature of an astronomical object 
moving through a medium. We examine the wake properties in detail and provide the obs ervable 
aspects in a quantitative way, by revisiting the previous semianalytic and numerical works (IKK07I ; 



Kim 



2010 ) and comparin g the features with the linear traje ctory counterparts discovered in analytic 



and numerical methods (|Ostrikerlll999l : iKim &; Kimll2009l ) 



The analytic formulae for the density wake in the previous work (reviewed in eq. [3] and [3]) 
do not give a clear physical insight of spiral- and/or arc-shapes. Hence it is difficult to deduce the 
properties of the wake relevant to the object properties. A further calculation in this paper leads to 
a purely analytic solution for the wake shape (eq. [5] ) , which allows detailed description of the wake 
shape. Simple geometrical approaches are adopted to provide the physical understanding of such 
shapes. Density enhancement profiles are also investigated. Since the empirical formulae based on 
the analytic analysis (eq. |ll]-|13j) outline the maximum and minimum density enhancement very 
successfully, these can be directly adopted for interpretation of future simulation and observational 
data. Our deeper understanding of the linear wake helps to understand nonlinear wakes shown in 
previous numerical simulations. Particularly, nonlinear aspects such as the increase of background 
density, the disappearance of inner arm boundary, and the appearance of a detached spiral shock 
are predictable from a close inspection of linear wake properties. 

Considering an object characterized by Bondi accretion radius rs and circularly orbiting at a 
distance r p with an orbital Mach number M p in a static uniform gaseous environment, our new 
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findings for linear wakes are as follows: 

1. In the orbital plane of the perturbing object, the linear wake shows excess density in the form 
of single-armed structure. For supersonic motion, it is confined by distinct boundaries of 
Archimedes spirals in the form of r/r p = tpA<i~ l +C with different constants C"s for the outer 
(always 1) and the inner (depending on M p ) boundaries. This expression for the explicit 
form of the spiral shape indicates the pitch angle of tan" 1 (M p r/r p )~ l and the arm spacing 
of 2irr p M p , from which we can deduce the orbital parameters, A4 P and r p . 

2. The vertical structure of the gravitational density wake is quantified for the first time. In 
meridional planes, the dense regions are bounded by two-centered concentric arcs forming 
crescent-like structures (M p > 1), partially overlapping each other for high orbital Mach 
number (M p > 4.6) cases. The angular size of the arcs are up to < 2tan -1 (.M2 — l) 1 ' 2 . 

3. The density enhancement a at the orbital center is always tbIt p . 

4. The density enhancement at the edges of the spiral arm pattern can be approximated to 
Oarm = «i + (4n + 2) a m i n with an integer number n depending on Ai p and some definitions 
of a\ (eq. [12] ) and a m ; n (eq. [IT]). The interarm density also increases to abater = (3ra— 2) a m ; n 
for supersonic cases (n > 1) and a m i n for subsonic cases (n = 0). Hence the density contrast 
between arm and interarm regions, (/9 a rrn — Pinter) /Pinter ~ aaim — center in the linear regime 
(aarm, center <C 1), is ~ ai + (?i + 4) u m \ n for supersonic cases and ~ u\ + a m in for subsonic 
cases. The numerical values for density contrast are approximately rs/r at large distance 
(r S> rv,), i.e., (0.7-2.9) rg/r for a wide range of orbital Mach numbers A4 P = 0-0.9 and 
1.5-10. 

5. The initially static fluid acquires velocity components in an approximate spherical form, 
starting from near the inner boundaries of high density regions and finally accumulating at 
the outer boundaries. 

6. The gravitational size of the object does not significantly affect the linear wake features insofar 
as the size is smaller than the interarm gap. 

As the accretion capability (rs/r p ) increases for a marginally supersonic orbital Mach number 
(A4 P ), the thermal pressure can be saturated at the junction of high density boundaries appearing 
as a tip in the linear analysis, so that the inner boundary extends beyond the outer boundary of 
the pattern. The induced wake reveals the following nonlinear behaviors: 

1. An appearance of a detached bow shock and a disappearance of the definite inner arm bound- 
ary in the highly nonlinear wake are mutually explained by an overpass of the inner boundary 
beyond the outer boundary. The shock standoff distance depends only on tb and M p , re- 
gardless of the adiabatic index 7. 
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2. Turbulent low density eddies are developed in the unstable process of producing the detached 
shock via the overpass of the inner boundary beyond the outer boundary. These eddies survive 
when the shock distance in the lateral direction is larger than rsM~ 2 , consistent to the linear 
trajectory counterpart in the nonlinear regime. 

3. The density enhancement a = rB/r p at the orbital center is preserved even in the extremely 
nonlinear regime, but the interarm density systematically increases proportional to this value, 
reducing the arm-inter arm density contrast. 

Using these theoretical diagnostics, we can derive the object mass, velocity, and distance from 
the orbital center with a given sound speed of the medium. In the case that concentric arcs are 
observed, the perturbing object is expected to be at the center of the circular arcs, potentially 
giving Tp. The angular size of the arcs corresponding to 2tan~ 1 (A^H — l) 1 ' 2 probes the orbital 
Mach number M p of the object of interest, or the corresponding velocity V p with the aid of an 
estimate of the sound speed c s . The object mass M p (or the accretion radius rs = GM p /c 2 ) can 
be roughly estimated from the density contrast ~ r^jr at a distance r from the center. High 
angular resolution observations are required to identify the broadening of arm edges (~ 2r s ) so 
as to estimate the object size r s . On the other hand, an observation of a spiral structure can 
provide information on r p and Ai p by measurement of arm spacing and arm width. An observation 
of smooth density distribution without sharp discontinuity provides an upper limit of the object 
speed (M p < 1). Finally, an observation of the arm structure with a missing inner boundary 
suggests the nonlinear wake condition, rs/r p (A4 p — l)" 1 > 0.1. 

From observations of spirals and/or arcs around stars, one may search for the possible existence 
of unseen companions, which are difficult to be directly detected either due to high obscuration 
by the dense environment or due to the bright glare of the central object. The presence of the 
companions may be inferred by their imprints on the environment at large distances from the 
star. For example, s piral structures have b een observed around AB Aur revealed in near-infrared 



2009) 



scattered emission (IFukagawa et al.ll2004l ). around AFGL 3068 in dust s cattered galactic light 



(IMauron fc Huggins! 120061 ) . and around CIT 6 in molecular line emission (|Dinh-V.-Trung fc Lim 



According to the result in this paper, the density contrast of a gravitational density wake due 
to a 10 Jupiter mass planet, orbiting at a distance of several AU, is expected to be of the order of 
10 % at a distance of a few hundreds AU from the cen tral star, where the sp iral arms are located in 



the circumstellar disk of the Herbig Ae star AB Aur (jFukagawa et al.l 12004 ). Although the density 
contrast of the order of 10 % is not significant, it may be distinguishable in this disk (rather than 
spherical) system. Whether such structures can be observed for low mass objects remains to be 
ascertained. And whether the density contrast due to the gravitational influence of planets has 
a similar order of magnitude in a disk system also needs to be checked. To clarify the claim of 



Lin et al.1 (J2006J) of forming giant planets as the origin of the observed spiral arms around AB Aur, 



it is necessary to carefully examine the density profile of the structure. 
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In the ca ses of AFGL 3068 and C I T 6, spiral structures are d etected on the scale of 1-10 
thousand AU (JMauron fe Hugginsl 120061 ; lDinh-V.-Trung fc Limll2009l ). which exclude the possibil- 
ity of planet wakes since the density contrast decreases with distance (< 1%). The large spiral 
structures in these evolved circumstellar envelopes a re likely to be caused by reflex motion of the 
mass losing star in com parable mass binary system (|Sokerlll994l ; lMastrodemos fc Morris!ll999l : lHe 



2007 



Edgar et al.ll2008l ). To date, AFGL 3068 is the only object that shows a clear spiral modula- 



tion in the circu mstellar envelope an d its companion star is indeed confirmed by a direct imaging 
in near-infrared (IMorris et al.ll200a ). The binary separation is about a hundred AU. Recently, 
Raga et al.l (|201ll ) estimated approximately 1 solar mass for the companion of AFGL 3068 from 
the spiral shape. With this companion mass, the gravitational wake is expected to contribute as 
much as 100 % of the envelope density at an observing distance of the wake of a thousand AU, 
doubling the density in the arm structure. This density contrast of the companion wake may be 
comparable to the density enhancement due to the reflex motion of the primary star. Hence, in the 
central region (< 1" by adopting the distance of 1 kpc to AFGL 3068), which will be soon resolved 
by the Atacama Large Millimeter/submillimeter Array (ALMA), there could exist two separated 
arms caused by two different mechanisms, i.e., one structure due to the reflex motion of the mass 
losing star and one gravitational density wake of the companion star. A comparison between the 
two mechanisms is needed in a range of parameter space relevant to the object and background 
properties. 

We note that our analysis is limited to objects in circular motion in a static uniform isothermal 
gaseous medium. For comparison with detailed observations, our assumptions for the background 
medium and the object motion may need to be relaxed. However, the wake features described in 
this paper may suffice to provide a zeroth order approximation of the primary features expected 
from the gravitational wakes due to such objects. 
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ing the manuscript, and Oscar Morata and Woong-Tae Kim for encouraging and helpful comments. 
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Fig. 1. — Density enhancement a of the spiral-onion shell shaped density wake generated by 
a point mass object orbiting at a fixed distance r p with an orbital Mach number A4 P = 10 in an 
ini tially s tatic uniform gaseous medium after 5 orbital periods, calculated by a semianalytic method 



m 



KK07I . (a) In the orbital plane, the induced density wake possesses a condensed trailing arm 
bounded by outer (black dashed) and inner (black dotted) Archimedes spirals r/r p = ip/ ' M. p + C with 
C = 1 and —0.7, respectively, in the polar coordinates (r, ip) with ip in the clockwise direction tracing 
the spiral but opposed from the motion of the object. White circle represents the circular orbit of 
the object, currently located at the intersection of white dotted lines, (b) In a meridional plane, the 
high density area of the wake displays arcs being able to extend in angle of tan _1 (.Mp — 1) ' 2 from 
the orbital plane (black solid). White dotted line indicates the current location of the object as in 
(a), (c) The density enhancement profiles as a function of the distance from the center are displayed 
along + x (red), ~y (yellow), ~x (green), + y (blue), and z (cyan) axes. Notice that a = rs/r p at the 
center, where r# is the Bondi accretion radius of the object. The black solid line (a = a± + 14a m i n ) 
outlines the peaks at the spiral arm boundaries, the dashed lines (a = a m j n x 1,4,7,10 from 
bottom) show the change of minimum density enhancement in the orbital plane because of piling- 
up of delayed perturbations, and the dotted line represents well the density enhancement of a m \ n 
along the orbital axis. See text for more comprehensive explanation including the definition of a\ 
and a min . 



22 
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Fig. 2. — Same as Fig. QJi-b, but for the number of perturbations. The perturbing object has 



circled 5 times with an orbital Mach number of A4 P = 10. 
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Fig. 3.— Same as Fig. lb, but for (a) M p = 0.8 (b) M p = 2.0, (c) M p = 4.5, and (d) M p = 5.0. In 
each panel, the n + 1 dashed and one solid lines represent a = a m in x 1, 4, • • • , 1 + 3n from bottom 
and a = ot\ + (4n + 2) a m i n , respectively, in the orbital plane with n of (a) 0, (b) 1, (c) 1, and (d) 
2. See text for the dependence of n on the orbital Mach number Ai p , and the definitions of a m i n 
and a\. 
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Fig. 4. — Density enhancement a of the spiral-onion shell density wake driven by an extended 
perturber with the gravitational softening radius of r s = 0.01r p in the unit of the orbital radius 
r p , calculated by a numerical simulation using an adaptive mesh refinement FLASH code with the 
maximum refinement level 8 corresponding to the spatial resolution of 0.001r p (as indicated by thin 
black boxes including 64 x 64 x 32 meshes inside). See Fig. [T] for the details. The thick black boxes 
in (a) and (b) denote the regions plotted in Fig. [6J and the vertical dotted lines in (c) represent the 
gravitational size r s of the object. 
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Fig. 5. — Same as Fig. [H but for a more extended perturber with the gravitational softening radius 
of r s = 0.1r p , comparable to half (0.09r p ) of the interarm spacing in the case of a point mass 
perturber. 
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Fig. 6. — Velocity field (arrows) superposed on the map of density enhancement a (logarithmic gray 
scale) for a segment of high density structure in Fig. 0J observed (a) in the orbital plane (z = 0) and 
(b) in a meridional plane (y = 0). The maximum velocity in each panel approximately corresponds 
to the line segment at the bottom right corner with the magnitude of A4 = r^Alr 1 |r — r p | _1 x 
2 (r s /r p ) -1 ' 3 in Mach number, where rg, M. p , r p , and r s represent the Bondi accretion radius, 
orbital Mach number, orbital radius, and gravitational size of the perturbing object, respectively, 
and r ~ 2r p is adopted. 
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Fig. 7. — Vertical distributions of density enhancement a (logarithmic gray scale) and velocity 
field (arrows) for the wake created by a gravitational potential characterized by its Bondi accretion 
radius r# and softening radius r s = 0.1r p , where r p is the orbital radius of the gravitational object 
about origin in z = plane. The object moving with Mach number Ai p = 5 is currently located 
at (x,y,z)/r p = (1,0,0). (a) Under a weak gravitational potential (tb/t v = 0.01), the induced 
wake has the shape of crescents bounded by circular arcs of radii A-rMZ, 1 + r s , QttM~ 1 + r s , and 
SnAdp 1 + r s about the object position with broadening by ±r s so as to appear banded (dotted) 
and of radii 1ttA4Z 1 , 3tt A4~ , and 5ttA4 p 1 about the mirror point from the object with the same 
broadening (dashed). The material flows approximately from inner dotted lines to outer dashed 
lines, (b) A strong gravitational potential (rB/r p = 1.0) introduces nonlinear features of the induced 
wake, different from the corresponding linear wake in (a). For comparison, the inner dotted and 
outer dashed lines in (a) are copied. Both density and flow speed of nonlinear wake are monotonous 
compared to those of linear wake. The line segment at bottom left corner represents the fiducial 
magnitude of velocity vectors in Mach number, M. = rBA4~ 1 \r — r p \~ 1 x 2 (r s /r p ) -1 ' 3 with adopting 
r ~ 2r p . The spatial resolution is 0.04r p in this area, although the maximum resolution achieved 
around the object position is twice as high with 5 levels of refinement in FLASH code. 
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Fig. 8. — Density enhancement a of a nonlinear wake due to a gravitational potential, characterized 
by the Bondi accretion radius r# = 0.5r p and the softening radius r s = 0.05r p , in circular motion at 
distance r p with the Mach number A4 P = 2.2. The calculation is performed by FLASH code with 
the maximum refinement level 5 corresponding to a spatial resolution of O.Olrv,. (a) In the orbital 
plane, the nonlinear wake departs in shape from the corresponding linear wake of a single-armed 
spiral confined in the dashed and dotted curves, but its shock front follows the solid curve, which 
is the dotted line phase shifted by n. (b) In a meridional plane, the nonlinear wake exhibits the 
filled elliptical arcs below the solid line as the vertical extension limit of the linear counterpart, (c) 
The density enhancement profiles along the distance from the center in the orbital plane (colored 
lines) have the maxima outlined by the black solid line (a = a\ + 6a m i n ) similarly to the linear 
counterparts, but the minimum density enhancement follows the dotted line (a = rs{r + r p )~ l ) 
instead of the dashed lines (a = a m ; n and 4a m j n ) for the minimum profiles in linear regime. See 
text and Fig. [1] for more details. 



